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Abstract 

Systems in which a short-ranged attraction and long-ranged repulsion compete are intrinsically 
frustrated, leading their structure and dynamics to be dominated either by mesoscopic order or 
by metastable disorder. Here we report the latter case in a colloidal system with long-ranged elec- 
trostatic repulsions and short-ranged depletion attractions. We find a variety of states exhibiting 
slow non-diffusive dynamics: a gel, a glassy state of clusters, and a state reminiscent of a Wigner 
glass. Varying the interactions, we find a continuous crossover between the Wigner and cluster 
glassy states, and a sharp discontinuous transition between the Wigner glassy state and gel. This 
difference reflects the fact that dynamic arrest is driven by repulsion for the two glassy states and 
attraction in the case of the gel. 
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I. INTRODUCTION 



In general, attractive interactions promote ordered phases or condensates, whereas long 
ranged repulsions inhibit this tendency, fundamentally redefining the system's free energy 
leading to complex phases that break translational and/or rotational symmetry. Such meso- 
scopic ordering occurs in a very diverse range of materials , from the pasta phase in 
neutron stars |3|, highly correlated quantum Hall and strongly correlated electron systems 
such as high superconductors j^, 4| to classical systems jl| such as ferromagnetic films, 
diblock copolymers, colloids 5|, |6|, and biological systems. It has also been suggested 
that competing interactions can also cause frustration, leading to exotic nonergodic disor- 
dered states. In the above examples, the system is able to relax on mesoscopic lengthscales; 
it is rather rare to see metastable disorder at the local level. Such disordered states were 
however reported for Laponite suspensions, where electrostatic repulsions compete with van 
der Waals attractions, but the anisotropic particles and interactions make the situation 
rather complex 

Spheres with a hard core repulsion and an attraction have long provided a model which 
captures the essence of atoms and small molecules. Short-ranged attractions lead to gelation 
due to arrested phase separation 10|. At higher densities both hard-sphere and attractive 
glasses are found [Fig. 1( g)] llj. along with gels 12|. Long-range repulsions can lead to 



glasses at low densities 13|, and combined with short-ranged attractions, the behavior is 



very rich and complex. Indeed, many properties of biological materials may be connected, 
fundamentally, to a short-range attraction and longer-ranged repulsion due to electrostatic 
charging from immersion in an aqueous medium 15|. For example globular proteins are 
rather well described as spheres with short-ranged attractions and long-ranged repulsions 
[igI . While most biological systems are much more complex than spheres with competing 
interactions, it is clearly important to understand this seemingly simple addition to well- 
studied models of atoms, not least as it offers insight into transitions between metastable 
states. 

Since competing interactions lead to frustration between phase separation and homogene- 
ity, a characteristic lengthscale is often predicted from computer simulation, for example, 
periodic lamellae jo] or low-dimensional clusters of a specific size jl^l. These may then un- 
dergo hierarchical self-organization; in particular, clusters may themselves be implicated in 
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gelation [18| and undergo dynamical arrest to form a 'cluster glass' [19|. Despite recent de- 



velopments 
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2l| , experimental work in colloidal systems with competing interactions 



las so far found little evidence of periodic structures, although gels with novel structures 
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21| and low- dimensional clusters 20|, |22| have been seen. 



II. EXPERIMENTAL 



We consider spherical colloids (diameter a) immersed in a solvent, with a relatively 
strong, long-range electrostatic repulsion, and a short-range, tuneable attraction mediated 
by non-adsorbing polymer in which the strength of the attraction is set by the polymer 
concentration Cp, and the range by the polymer size, in this case the polymer-colloid size 
ratio g ~ 0.19 (see appendix). We determined the attractions in a very similar system 23| . 
and the magnitude of the electrostatic repulsions from fitting the structure of equilibrium 



fluids 



24| which gave a colloid charge number o 



to mode-coupling theory for Yukawa systems 



Z ~ 600 ± 200 (see appendix). According 
251 ] ■ for our parameters we expect a transition 
to a Wigner glass at comparable colloid volume fractions to those observed experimentally. 
Mixing the samples prior to imaging leads to a randomized initial state, after which the 
tuned interactions and colloid volume fraction (p lead the system to a (metastable) point on 
the state diagram (Fig. [1]). We sometimes see the formation of a Wigner crystal state at 
high (f) > 0.2 and low Cp; however, this is a rare event [see Fig. 1(f)] and usually the system 
forms a metastable glassy state before crystallization can occur (see appendix). Using real 
space structural and dynamical analysis, we find a state diagram dominated by three types 
of glassy states with different disordered structures: Wigner glassy state, cluster glassy 
state, and gel. We study the transitions between these states, and reveal the nature of these 
transitions and its link to the interactions leading to slow dynamics. 



III. RESULTS AND DISCUSSION 



We begin by presenting the state diagram in Fig. [T], which underlines the extent to 
which the system is dominated by dynamically arrested states. A low-density colloidal fluid 
(0 = 0.002) where the system is ergodic is shown in Fig. W^c). Increasing the volume 
fraction to (p = 0.016 results in a glassy state where the slow dynamics is driven by the 



3 







Figure 1: State diagram of the investigated system and 2D and 3D structures, (a) cluster glassy 
state, 4> = 0.04, Cp = 8.39 g/1. (b) Wigner glassy state, </> = 0.047, Cp = 0. (c) fluid, 4> = 0.002, 
Cp = 0. (d) gel, 4> = 0.152, Cp = 4.10 g/1. (e) Wigner glass, (/) = 0.372, Cp = 0. (f) crystal, 
4> = 0.372, Cp = 0. Bars are 20 /xm. Boundaries are guides to the eyes. Increasing grey shading 
represents increasing amount of clusters. Thick arrows denote paths 1-3 described in the text. 
Cluster glassy state is denoted by squares, gel by circles, crystal/glassy state by hexagons and 
fluid by triangles. Grey arrows indicate states for the images (a)-(f). (g) State diagram of a 



111 ], showing two glassy states at high colloid 



system without long-ranged repulsive interactions 
concentration, (h)-(j): 3D structures of (h) Wigner glassy state, (i) cluster glassy state, and (j) 
gel. White particles have no neighbors, otherwise colors denote connected regions. 
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long-ranged electrostatic repulsions [Figs. [U^b), (e), and (h)]. We thus term this state a 
'Wigner glassy state'. At low and higher Cp, we see the formation of clusters and term this 
state a 'cluster glassy state' [Figs. [Tt^a) and (i)]. Meanwhile, increasing both Cp and (p results 
in a gel which we define through percolation [Figs. [U^d) and (j)], and appears dynamically 
arrested [Fig. Et^a)]. A typical state diagram for a colloidal system without electrostatic 
repulsions is schematically shown in Fig. 1(g). We see striking differences in the structure 
of the state diagram between systems with and without electrostatic repulsions (see also 
appendix). With electrostatics the ergodic region is very narrow, around ~ 0, and the 
glassy states dominate the state diagram. There are three glassy states with very different 
structures: Wigner and cluster glassy states and a gel state. We also found three types of 
transitions between these states accompanied by local structural changes. 

Before discussing the structure of these three states in more detail, let us consider the 
dynamics. Mean squared displacement (MSD) measurements, in which the colloids are 
tracked in two dimensions (2D), are shown in Fig. |2j What is clear is that the (ergodic) 
fluid we see at low appears to exhibit diffusive behavior, and the gel appears dynamically 
arrested, within the accuracy to which we can track the particles (100 nm or (t/20). The 
other states show extensive non-diffusive behavior, yet they do not reach a clear plateau on 
the experimental timescale although we track the particles for up to 20 hours for the cluster 
glassy state and up to 2 hours for the Wigner glassy state. We are limited in particle tracking 
at long times due to particle loss (particles diffusing out of our volume), bleaching and drift, 
along with residual sedimentation. We estimate the cage size as the width of the first and 
second peaks in g{r) for the Wigner and cluster glassy states respectively [see Fig. 3(a)]. In 
the case of the gel the cage size is taken as the bond length, O.lOa. Concerning the cluster 
glassy state, similar sub-diffusive behavior has recently been seen in computer simulation 



191]. We note that such behavior is expected from the fact that clusters can rotate as rigid 



bodies [see Fig. 2(b)]. This strong decoupling between translational and rotational motion 
may be unique to a cluster glassy state. In the case of the Wigner glassy state, based upon 
simulation work in a similar system 26|], we believe that our experimental timescales limit 
full access to relaxation phenomena. In any case, the cage size is not reached [see Fig. 2(a)]. 
The soft nature of the interactions is another factor which prolongs the sub-diffusive regime 
and thus makes a rigorous confirmation of nonergodicity very challenging. The non-diffusive 
behavior is suggestive of the glassiness of these states. Since we do not find a clear plateau 
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Figure 2: (a) Mean square displacements. Crosses=fluid ((/> = 0.002, Cp = 0), green trian- 
gles=cluster glassy state {(p = 0.051, Cp = 5.16 g/1), red circles=gel {(p = 0.275, Cp = 4.88 g/1) 
and blue squares=glassy state {(j) = 0.069, Cp = 0). r=characteristic diffusion time (see appendix). 
Dashed lines corresponding to the cage size are plotted for gel (0.19(t), Wigner state (1.9cr) and 
cluster state (2.1cr), respectively, (b) Confocal microscope images of a cluster spinning in it's 
'cage'). Time unit=r, image width=19.5 /.im. 

experimentally, we term these 'cluster glassy state' and 'Wigner glassy state', respectively, 
although the states may be regarded as nonergodic at least practically. 

The state diagram dominated by these three glassy states illustrated in Fig. [1] yields 
three transitions. It is a commonly held view that dynamical arrest is accompanied by 
little structural change. While literature on transitions between glassy states is scarce, the 
structural change in the attractive-repulsive glass transition at high packing fraction is rather 
subtle j^^]. Conversely, all states identified here are characterized by their structures [see 
Figs. 1(h), (i), and (j)]. Our analysis shows a considerable variety of metastable structures 
in the response to small changes in parameters. 
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Figure 3: Dependences oi g{r) and (m^) on (j) and Cp. (a) (path 1) shows the change in g{r) across 
the Wigner-cluster glassy state transition. Vertical solid line indicates the peak at contact (r =o"), 
and dashed line indicates the first peak of g{r) of the Wigner glassy state, corresponding to the 
average interparticle distance, and the second peak of g{r) of the cluster glassy state, corresponding 
to the first shell of clusters, (b) (path 2) from a Wigner glassy to a gel state, and (c) (path 3) from 
a cluster glassy to a gel state. Cp is given in units of g/1. Data offset for clarity, (d) The change 
in (?7T,^) along path 1 (at cj) = 0.051), (e) along path 2 (at cp = 0.151), and along (f) path 3 (at 
Cp = 7.58 g/1). Dashed lines and shading are guides to the eye. 

Let us now enquire as to the nature of the transitions between these glassy states. The 
transition between Wigner and chister glassy states (path 1 in Fig. [1]) is considered in more 
detail in Figs, and (d), which show the radial distribution function g{r) and the variance 
in the cluster size distribution (m^), where m is the number of particles per cluster. Rather 
than a sharp transition, the Wigner glassy state is unaffected by weak attractions; until the 
polymer concentration exceeds around Cp = 4 g/1 there is no response in the size distribution. 
At higher Cp, there is an increase in (cluster) size, yet our data suggest that it occurs rather 
gradually, i.e., passing from the Wigner to the cluster glassy state is a crossover rather 
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than a sharp transition: particles start to form small clusters above a certain threshold 
Cp and their size gradually increases with an increase in Cp. The gradual development of 
the peak at contact (at r = a) in the g{r) [Fig. ^a)] further supports this observation. 
Note that the emergence of the peak at contact is a direct consequence of attractions. In 
equilibrium, the repulsive Wigner state (comprised of monomers) is expected to be a crystal. 



The same can apply to states of monodisperse clusters [28[, however here the clusters are 
very polydisperse [see Fig. dl^a) and appendix], which leads to self-generated disorder, and 
intrinsically suppresses crystallization. This source of disorder is a consequence of a two-level 
organizational hierarchy, of colloids forming clusters and then clusters forming a glass which 
may be characteristic of a system of competing interactions with different lengthscales. 

Proceeding to the transition between the Wigner glassy and gel states (path 2 in Fig. [1]), 
we find a rather different scenario. Raising the colloid volume fraction to </> = 0.15, the g{r) 
[Fig. 3(b)] again shows the development of a peak at contact (r = a) around Cp = 4 g/1. We 
recall that at lower 0, at a similar polymer concentration, clusters began to form [Fig. |3]^a)]. 
For (f) = 0.15 this yields percolation [Fig. [U^j)], and a sharp transition to a gel state [Fig. 
3(e)], accompanying a strong increase in (m^) by about three orders of magnitude. This 
is markedly different from the gradual continuous increase in (m^) for the Wigner-cluster 
glassy states [see Fig. 3(d)]. 

What happens in the case of the transition between cluster glassy and gel states? Path 
3 in Fig. [T]is shown in Figs. [3t^c) and (f). Unlike the previous transitions, here the polymer 
concentration is fixed at Cp = 7.4 g/1. In fact rather little change in the local structure is seen 
in the g{r) analysis [Fig. Et^c)], and the variance in the cluster size increases continuously. 
Our results indicate that, in contrast to paths 1 and 2 in Fig. [T]the cluster glass-gel boundary 
is delineated by a percolation transition, rather than by local structural changes. 

We also observed a novel aging mechanism in the cluster state. While aging is usually 
thought of in dynamic terms, here we present a structural mechanism for aging. This 
concerns the transition from a state of high potential energy, a tightly bound cluster of 
charged particles, to a smaller cluster, via cluster fission. We never observed any cluster 
fusion, and particle tracking shows a continuous rise in the population of small clusters as 
a function of time. The emission process occurs in less than 1/100 of the characteristic 
diffusion time, a much faster timescale than structural relaxation even in the absence of 
slow dynamics. This phenomenon is illustrated in Figs. Hl^a) and (b). This cluster state 
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Figure 4: Aging mechanism of a cluster glassy state, (a) An emission process from an m = 5 to 
m = 4 cluster, as shown by arrows. Time t is expressed in units of r/1000 and the image width is 
14 /im. (b) Separation of the emitted particle and cluster, as defined by arrows in (a), as a function 
of t. 



suffers from a complex minimization problem involving the spatial distributions of colloids, 
counterions and polymers. The coupling between these three variables makes the potential 
energy landscape rather complex. This fission process allows us to directly observe a kinetic 
path from one local minimum to a neighboring minimum with a lower energy. Details of the 
cluster distribution and its time-evolution are available in the appendix. 

Some comments on the colloidal interactions are in order. We have determined the colloid 
charge at low </> and Cp = 0. Using the parameters we calculate for the colloid charge, for this 
system we expect a potential at contact of several hundred k-oT, substantially in excess of the 
attractive forces mediated by the depletion attraction, ~ 30 k-oT for Cp = 4 gl~^ at which 
we find clustering and gelation (see appendix). We believe that the charge is associated 
with the free surfaces of the colloids, and that clusters may have a markedly lower charge 



per colloid than do the isolated particles 



29j. Calculations show that for such anisotropic 



charge distributions, the electrostatic repulsion is much reduced, enabling the clustering and 
gelation we observe jsO ]. 



IV. CONCLUSIONS 



In closing, we found that the introduction of relatively strong, long-ranged repulsions 
to a colloidal dispersion with short-range attractions, generates novel glassy states, such 
as the cluster glassy state, and drastically transforms behavior, at that most fundamental 
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of levels: the ability of the system to relax locally. For a system of short-ranged attractive 
interactions, without strong, long-ranged repulsions, it is possible to distinguish the repulsive 
hard-sphere glass and the attractive glass, as depicted in Fig. [U^g) ll|, but there the change 



in the structure is rather subtle 



27|. In contrast, the transitions between cluster glassy and 



gel states and between Wigner glassy and gel states both accompany significant structural 
changes from single particle to compact clusters and a percolating network, respectively. The 
observed behavior can be interpretted as the interplay of a Wigner glassy state dominated 
w the long-range repulsions with formation of clusters driven by short-range attractions 



17l | and their percolation for gelation. Like the neutral system, there appear to be two 
drivers of arrest: the gel is driven by attractions and both Wigner and cluster glassy states 
are driven by repulsion. The effect of long-ranged repulsions on the state diagram of a 
system with short-ranged attractions is an interesting fundamental problem; at intermediate 
repulsions a cluster fluid emerges 



20 



22|. Since protein solutions are also understood to 



exhibit comparable interactions, we conclude that it is at least reasonable to suppose that 



equihbrium cluster phases may be found 



16 



3l|, and perhaps even cluster glasses. 
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Appendix 

Samples and experimental methods. We used poly(methyl methacrylate) (PMMA) 
colloids sterically stabilized with polyhydroxyl steric acid. The colloids were labelled with 
the fluorescent dye rhodamine isothiocyanate and had a diameter a = 1.95 fim with around 
5% polydispersity. Solvents and polymers were used as received. The polymer used was 
polystyrene, with a molecular weight of M„ = 2.06 x 10^. To closely match the colloid density 
and refractive index we used a solvent mixture of cis-decalin and cyclohexyl bromide. Due 
to the refractive index matching, the van der Waals interactions are reduced to a fraction 
of the thermal energy /cbT {ks'- Boltzmann's constant; T: temperature) and neglected. We 
confined the samples to thin (100 /im) capillaries which provided access to the entire sample. 
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but confirmed, with experiments using larger 500 fim cells, that this confinement has little 
effect upon the behavior. Prior to each experiment, sample vials were washed with copious 
quantities of ethanol, and carefully dried. Dry colloids were then added, and dispersed in 
solvent, with the polymer solution added last. The data was collected on a Leica SP5 confocal 
microscope, fitted with a resonant scanner. We determined the co-ordinates of each particle 



with a precision of around 100 nm (23|. Connectivity was determined by setting the bond 
length equal to the sum of the range of the attractive interactions and the tracking error. 
Our conclusions are insensitive to the exact value of the bond length. The time for a colloid 
to diffuse one radius is given by the Stokes-Einstein relation as r = Srcqa^ / (Ak-BT) where 
rj is the viscosity. This we measure independently as a function of polymer concentration 
using a Rheologica Instruments Stress Tech rheometer. 

Determining the interactions between the colloids. The behavior we have observed is 
driven by competing attractive and repulsive interactions. We shall begin our discussion of 
the potentials by considering the attractive depletion interaction. We have used a polymer 
of Mw = 2.06 X 10^. Here we used a 'good' solvent, and in detailed studies using the same 
solvent-polymer mixture, we found good agreement with the Asakura-Oosawa (AO) theory 



32|, assuming polymer swelling such that the radius of gyration Rq is increased by 35%. 
We therefore assumed a similar degree of swelling here, such that Rq ^ 190 nm, leading 
to a polymer-colloid size ratio q = 2Rq/(J ~ 0.19, which sets the range of the attractive 



interaction. For the dilute polymer regime, little deviation is expected from AO theory 32| . 



even in the case of some polymer swelling (non-ideality) 33|. Furthermore, although we add 



sufficient polymer that some of our samples are just into the semi-dilute regime, relatively 

. At the level of this 



little deviation is expected from AO theory for these parameters 
work, therefore, we treat the attractions with AO theory which leads to a pair interaction 
between two hard colloidal spheres in a solution of ideal polymers which reads 

oo for r < a 



/3uAo{r) 



6 

3r 



(1) 



><{l-2{if^ + 2(Tw} ioTr>a<a + {2Ro) 
for r > a + {2Rg) 

where (3 is l/ksT. The polymer fugacity zp^i is equal to the number density ppr of ideal 
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polymers in a reservoir at the same chemical potential as the colloid-polymer mixture. This 
corresponds to a contact potential of 7.0 k^T per g/1 of polymer reservoir concentration. 
Determining the electrostatic interactions in non-aqueous colloidal dispersions is a non- 



trivial task 



23 



34l ]. as the charge and ionic strength are very small, although they still 



lead to significant interactions. The interactions have nevertheless been found to be well 
described by the screened Coulomb or Yukawa interaction, which reads 



f3uYUK{r) 



oo 



/3e 



exp(— K(r— (t)) 
r/a 



for r < (T 
for r > a 



(2) 



where r is the center to center separation of the two colloids. The contact potential is given 
by 

Pe= ^^ ,^\J -^, (3) 



(1 + Ka/2)2 a ' 

where Z is the colloid charge, k is the inverse Debye screening length and Zb is the Bjerrum 
length. 

One may thus assume interactions between the colloids take the form of a short-range 
'depletion' attraction and long-ranged screened Coulomb repulsion, and a nearly hard core. 
Both have been measured in nearly identical systems and found to be in good agreement with 



theory 
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24| . Here we have determined the colloid charge and ionic strength by comparison 



of the measured g{r) with one obtained from Monte Carlo simulation according to a Yukawa 



potential 



24| in the ergodic fluid part of the state diagram 



4.8xl0-4±2xl0-^ c 



0) 



tohe Z ^ 600 ±200 e where e is the elementary unit of charge, and the inverse Debye length 
K,a ~ 1±0.4. The radial distribution function g{r) fitting is shown in SFig. 1. This method is 



found to be consistent with electrophoretic measurements 



24l ]. The values we obtain for the 



colloid charge are somewhat (up to 25%) higher than those quoted previously in literature 



3, 



34| 



or these parameters mode-coupling 



25| . However, this analysis poses a 



where samples were prepared using dry colloids 
theory predicts vitrification at higher volume fractions 
significant question: the contact potential between two such colloids is (3e ~ 1000. This is 
much greater than the strength of the attraction induced by the polymer. Throughout the 
measurements we find that clustering of gelation occurs at a polymer concentration around 
Cp = 4.0 gl~^, which corresponds to a contact potential of around 30 [Eq. ([1])]. 

We are reasonably confident about our understanding of the attraction induced by the 



polymer 



23| . We speculate that the charge is in fact different between the isolated colloids 



12 



upon which measurements of the charge Z are based and colloids found in clusters, consistent 



with behavior we observed in a similar cluster-forming system 
the literature 
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20 



29| and values quoted in 



2l|. Alternatively, non-pairwise additivity with these long-ranged 



repulsions leading to a reduction in apparent repulsion may play an important role 35 |. 
In all of these cases, the quoted values for the repulsive interactions are stronger than the 
depletion induced attractions. 

Pointers to an understanding of the underlying cause for this discrepancy may be found 
in measurements that indicate a reduction in charge per colloid upon clustering 29j . More- 
over, if this charge were inhomogeneously distributed, a further reduction in repulsions is to 



be expected 



sol ] upon which clustering due to polymer- induced attractions is not unreason- 



able. Further work is in progress to elucidate the charging mechanism, and possible charge 
anisotropy, however, we note that the weak levels of electrostatic charging in these apolar 
systems severely hamper a full understanding. 

However, the strength of the depletion interaction, as gauged by its strength at contact 
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36|. 



[Eq. ([T])] is rather stronger here than in systems where the charge is screened 
Comparing the minimum polymer concentration required to form a gel, we find around a 
contact potential of fiuAo ~ H- When the electrostatic repulsions are screened one finds 



/3m AO ~ 3 for attractions o: 
shorter-ranged attractions 



comparable range (q=0.18) |36|, rising slightly to /3mao ~ 5 for 



10| . Previous work on charged systems with a much short-ranged 



attraction [ q = 0.04) suggests a threefold increase in the contact potential /3uao required 
for gelation 2l|] relative to systems where the charge is screened. 
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SFig. 1: Analysis of g(r). Radial distribution function. Dashed red line corresponds 
to a colloid charge Z = 400, solid cyan to Z = 800. We assumed that the Debye screening 
length was dominated by the colloidal counterions, in other words the system is close to the 
salt-free limit. This leads to a fitting which depends solely upon Z. Lower values of Z gave 
poor fits, higher values of Z led to crystallization. 
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SFig. 2: g(r) for repulsive glass. Decay of the spatial density correlation at long range 
suggests an absence of crystallization and that the state is a repulsive Wigner glass. Polymer 
concentration Cp is given in units of g/1. 



On crystallization. Occasionally we found that, rather than a Wigner glassy state [Fig. 
1(e)], the system in fact crystallized [Fig. 1(f)]. That similar samples may either crystallize 
or vitrify su gge sts multiple paths through the free-energy landscape, as found in some clay 



suspensions 



37| . Alternatively, the systems size we study allow the possibility of a variation 



in nucleation time between different experiments, leading to the observation of crystallization 
in a limited set of samples. 

The fact that the crystal forms at all also suggests that the increase in (j) may promote 
ordering and reduce the glass-forming ability of the system. We never observed an ordered 
phase for (p < 0.2, this is further supported by the radial distribution function in SFig. 2 
which shows an absence of long-ranged order for the repulsive glass. That we only found a 
Wigner glass at low 0, an apparently shallower quench, but occasionally a crystal at higher 
0, is interesting, and may be explained by a stronger screening of the electrostatic repulsion 
due to a higher concentration of counter ions at higher 0, leading to 'harder' interactions 
which promote crystallization. 
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m t [h] 

SFig. 3: Cluster distribution and aging mechanism, a, shows the variation in cluster 
size distribution over time. The inset shows the distribution for clusters of size m < 5. n 
denotes the number of clusters, the number of coordinates (particles) found, m gives the 
size of clusters, b, shows the evolution of different cluster sizes, where An(t) = n(t) — n{0). 
Note the increase in the population of m = 6 at early times. 

On the absence of mesoscopic ordered structures. Here we note that a consider- 
able literature exists, predicting pattern-formation and lamellae-like phases in systems with 
competing interactions 0, We see little periodicity in the system studied here. In the 
case of the cluster glass, this is likely due to extremely slow equilibration to a monodisperse 
cluster phase which could then potentially form a cluster crystal 28|. Our experiments run 
for up to two weeks, and there was no indication of crystalline ordering, moreover the cluster 
distribution was still very polydisperse on this timescale. The behavior up to 18 hours is 
shown in SFig. 3. 

In the gel, on the other hand, we believe that the bicontinuous network is slow to break 



and re-arrange due to the connectivity [38j, as would be required to form a lamellar type 
structure. Another important difference in comparing the gel state with for example, mi- 
croemulsions, lies in the arrested nature of the 'colloid'-rich phase. Microemulsions are thus 
able to reorganize locally in a way that a system with short-ranged attractions, is not. 



Dependence of the state diagram of multiple dynamic arrested states upon the 
strength of electrostatic interactions. The state diagram of multiple dynamic arrested 
states for varying strengths of the electrostatic interactions is shown in SFig. 4. Upon 
reducing the strength of the electrostatic interaction, one expects to recover the behavior 
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of hard spheres with a short-ranged attractive interaction ll|. That is to say, fluid, gel, 
attractive and hard sphere glasses. Here we discuss the differences this phase diagram 
relevant to neutral colloids to the strongly charged case presented here. Note that we 
consider the case of no added salt such that the electrostatics are unscreened. 

The effect of increasing electrostatic repulsions is as follows. When the colloid charge Z is 
small, a fluid (blue), gel, attractive glass (yellow) and repulsive glass (red) states are found 



Our study suggests that increasing the electrostatic repulsions leads to the system 



studied here, with a significant extension of the repulsive glass to low densities, such that we 
find a direct repulsive glass-gel transition. This is quite different from the fiuid-gel transition 
in the case of small Z. Moreover the enhancement of repulsions due to the electrostatics 
leads to an increase in the amount of polymer required for gelation. 

We furthermore see the emergence of a cluster glassy state (green), which is not found 
for neutral systems. Nonergodicity stems from connectivity in the case of the gel and caging 
for single particles and clusters, respectively, for the repulsive and cluster glasses. For a 
high colloid volume fraction, there may also be an attractive glass state, although such a 
high volume fraction region was not investigated in this study. We emphasize that for large 
Z an ergodic state exists only for a very low colloid volume fraction and most state points 
form non-ergodic disordered states. It is an intriguing, fundamental question how the state 
diagram evolves between these two extremes. 
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SFig. 4: 3D state diagram. The state diagrams for two extremes, negligible and strong 
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electrostatic interactions {Z = and large Z) are indicated. The former has been established 
in 111], while the latter is a schematic of Fig. 1. See the text for details. 
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